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A new method for determining physical parameters of 
fundamental mode RR Lyrae stars from multicolour light 
curves 
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ABSTRACT 

We present a new method for determining physical parameters of RRab variables ex- 
clusively from multicolour light curves. Our method is an Inverse Photometric Baade- 
Wesselink analysis which, using a nonlinear least squares algorithm, searches for the 
effective temperature (T c ff) and pulsational velocity (V p ) curves and other physical 
parameters that best fit the observed light curves, utilising synthetic colours and bolo- 
metric corrections from static atmosphere models. The T e g and V p curves are initially 
derived from empirical relations then they are varied by the fitting algorithm. The 
method yields the variations and the absolute values of the radius, the effective tem- 
perature, the visual brightness, and the luminosity of individual objects. Distance and 
mass are also determined. The method is tested on 9 RRab stars subjected to Baade- 
Wcsselink analyses earlier by several authors. The physical parameters derived by our 
method using only the light curve data of these stars arc well within their possible 
ranges defined by direct Baade-Wesselink and other techniques. A new empirical re- 
lation between the Iq magnitude and the pulsational velocity is also presented, which 
allows to construct the V p curve of an RRab star purely from photometric observations 
to an accuracy of about 3.5km/s. 

Key words: stars: horizontal branch - stars: variables: other - stars: oscillations 
(including pulsations) - methods: data analysis - techniques: photometric - 



1 INTRODUCTION 



The light variations of pulsating variable stars contain 
much useful information both on the physics of the bright- 
ness changes and also on the global parameters of the 
objects. Concerning RR Lyrae stars, this information can 
be retrieved using e.g., dif ferent varian t s of the Baade- 
Wesselink (B-W ) ana l ysis (iBaadel Il926t IWesselinkl 1 19691 : 
iManduca fcBeH Il98ll : iGautschvl 1 19871 : iFernlevI Il990h or 



empirical relations between the Fourier parameters of the 
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). However, there are cases, for example, 


Blazhko 



RR Lyrae stars showing light curve variation when we are 
encountering problems with both methods. The empirical 
light curve - physical parameter relations are proved to be 
valid only for RR Lyrae stars with reasonably stable light 
curves i.e., for unmodulated or weakly modulated stars. A 
good example for the latter is RR Gem, which is one of the 
calibrator objects of these relations, wh ile its light curve is 
weakly but unambiguously modulated (| Jurcsik et alj [20051 : 



ISodor et al.ll2007T ) . In spite of its small amplitude light-curve 
modulation, RR Gem fits well all the empirical light curve 
- physical parameter relations. The Fourier parameters of 
the mean pulsation light curve of a large modulation am- 
plitude Blazhko star, however, differ significantly from the 
parameters calculated for the li ght curves at any pha se of the 
modulation (see e.g. MW Lyr, Jurcsik et al.ll2008al ). which 
makes the inference of the empirical relations ambiguous for 
these Blazhko variables. Applying the empirical formulae on 
the Fourier parameters of MW Ly r in different phases of the 
modulation (data from Table 7 in I Jurcsik et alj[2008af ) vari- 
ations of 0.45 dex in [Fe/H] are obtained during the Blazhko 
cycle. It would be difficult to find a physical explanation for 
such a metallicity variation. Moreover, the magnitude aver- 
aged absolute brightness changes calculated from the Fourier 
parameters are 0.025 mag larger than what is in fact directly 
observed. Therefore, at present the light curve-physical pa- 
rameter formulae are not appropriate to follow the changes 
in the mean physical parameters during the Blazhko cycle. 

The B-W analysis is capable for determining physical 
parameters of Blazhko stars for any phase of the modulation 
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if simultaneous V la d and multicolour photometric observa- 
tions were available. Thus, any changes in the mean physical 
parameters could be, in principle, detected. However, a prac- 
tical problem is that for obtaining such results simultaneous 
extended spectroscopic radial velocity and multicolour pho- 
tometric observations were needed with good phase coverage 
of both the pulsation and the modulation cycles. For such 
an investigation at least one order of magnitude more Kad 
data points were needed than are available today. 

Up to now, there are only about two dozens of unmod- 
ulated field RRab stars, about one dozen RRab stars in 
globular clusters, and no Blazhko variable which have si- 
multaneous photometric and V la ,d observations suitable for 
B-W analysis. On the contrary, extended and accurate mul- 
ticolour observations of many modulated and unmodulated 
RRab stars have become available from our CCD observa- 
tions at the Konkoly Observatory in the past several years 
jSodorl 120071 ). To extract physical information exclusively 
from multicolour light curves of RRab stars we have devel- 
oped a method that fits the observed magnitudes and colours 
with appropriate physical parameters and their variations. 
The main advantage of the method is that all these are de- 
rived exclusively from photometric observations without any 
need of spectroscopic measurements. 

In this paper our method is introduced and is tested on 
9 non-modulated RRab stars. The application of the method 
on a large modulatio n amplitude R R Lyrae star, MW Ly r 
ijjurcsik et alJl2008eJ ) . is presented in ljurcsik et al.1 (|2008bh . 



2 THE METHOD 

The basic idea of the Inverse Photometric Baade-Wesselink 
method (IP method) is to look for the pulsational velocity 
(Vp(ip)) and effective temperature (T e ff(y>)) curves (<p de- 
notes the pulsation phase), the distance, and the mean ra- 
dius of the star which best reproduce the input light curves 
utilising static atmosphere model results. 

The inputs of the method are the Fourier series fitted 
to the multicolour BVIc light curves and the metallicity 
([Fe/H]). The input Fourier curves may represent either the 
light curves of an unmodulated RRab star, the mean pul- 
sation light curves of a modulated star, or the light curves 
of a modulated RR Lyrae star at a certain Blazhko phase. 
To obtain correct results the extinction corrected apparent 
visual brightness and the intrinsic colours have to be known. 
However, as it is shown in Sect. [4721 the latter requirement 
is not always crucial, since the colour zero points can be de- 
termined by the method itself if good quality light curves 
are available. 

For using model atmosphere tables the value of [Fe/H] 
and the instantaneous values of T e g(tp) and the effective 
gravity (\ogg e g((p)) are needed. The latter is determined 
from the V p (ip) curve, the mass (SOT) , and the average ra- 
dius (-Ro) of the star, as follows. The instantaneous radius 
(R(tp)) and the radial acceleration (a(v)) are computed from 
V p (<p) as the first integral and the first derivative, respec- 
tively. Knowing these variables, g c s[}p) are then calculated 
according to the formula: 

ges(<p) = Gm/R(ip) 2 + a(ip) (1) 



where G is the gravitational constant. The mass satisfies the 
modified form of the pulsation equation: 

log SOT = -2.237+2.56(logK)o-1.5381ogP-0.054(logT eff )o(2) 

which has been derived from Eq. la of iMarconi et al.l 
(|2003T ) by substituting the mean logarithm of the luminosity, 
(log L)o with: 

(log L) = 4(log T cff )o + 2(log R) - 15.048. (3) 

Note that SOT, L, and R are in Solar units and '0' index 
denotes pulsation-cycle-averaged values. P is the pulsation 
period of the variable in days. 

Synthetic colours, (B — V)(ip) and (V— Ic)(<p), and bolo- 
metric correction, BC{ip), corresponding to [Fe/H], T e ff(ip) 
and log fl c ff(<p) are taken from the static atmosphere model 
tables of ICastelli fc Kuru cz (2003). We applied linear inter- 
polation in [Fe/H] to obtain a resolution of 0.05 dex between 
the 0.5 dex steps of the original tables. The original resolu- 
tion of the and logg grid are 250 K and 0.5 dex, respec- 
tively, which have been interpolated by spline functions. The 
extinction corrected apparent magnitudes of the star is then 
derived from: 

V(tp) = -10logT(<p)-5\ogR(<p)-BC(<p)+51ogd+37.36(4) 

where d is the distance in pc. 

An implementation o f the Levenberg -Marquardt non- 
linear fitting algorithm by iLourakij (|2004h is applied to fit 
the derived light and colour curves to the observations with 
the alterations of the initial T e s(<p) and V p (ip) curves, d, 
and Ro. To obtain physically plausible results continuity, 
smoothness, and periodicity of the variations of the physical 
parameters have to be maintained during the fitting process. 
In order to attain these, the initial T e ff(<p) and V p {ip) func- 
tions are defined as smooth and continuous funcions which 
are then modified by the fitting algorithm with the addition 
of appropriate order Fourier series. 

RR Lyr atmospheres differ from static atmospheres in 
many respects. Just to mention some of the most problem- 
atic issues, in the dynamic atmosphere of large amplitude 
pulsating variables shock waves propagate, local thermody- 
namic equilibrium does not fulfil (NLTE), the atmosphere 
is not thin compared to the radius of the star while plan- 
parallel models are calculated, the turbulent velocity (v t ) 
also varies during the pulsation cycle, etc. However, as at- 
mosphere models taking into account all these complications 
(or at least some of them) have not been calculated yet for a 
wide enough range of parameters, all the methods which are 
applied today to define the physical parameters of pulsat- 
ing variables from photometric and spectroscopic data relies 
on static atmosphere models to some extent. The applica- 
bility of static atmosphere models for dynamic atmospheres 
of pulsating RR L yrae stars has already been discussed by 
many authors (e.g. I Liu fc Janesl Il99d : I Jones et al.lll992T l. It 
has been shown that the vicinity of the minimum radius 
phase is that part of the light curve where the validity of 
static models is the most questionable. Notwithstanding, in 
the IP method the whole phase range of the pulsation is 
used in order to maintain the continuity of the curves. The 
smoothness and continuity of the curves are especially im- 
portant because, due to the dynamic nature of the atmo- 
sphere, and to possible differences of the temperature scales 
the different colours corresponding the input light curves are 
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sometimes out of the ranges of the static atmosphere grid 
e.g., on the (B — V)— (V — Ic) plane which would otherwise 
introduce discontinuity in the output curves. Since we do not 
use the input colour and magnitude values at any phase di- 
rectly for calculating physical quantities, problems that arise 
from inappropriateness of the static atmosphere models at 
certain phases does not affect our results directly. Our in- 
verse method ensures that the solution curves are completely 
within the ranges of the model grid. Where the static model 
is less appropriate a poorer fit is expected. 

2.1 The initial T c k(<p) curve 

The (V — Ic) colour is a good temperature indicator both 
theoretically and e mpirically. We use the polynomial trans- 
formation given bv lBessel et al.l (|l998l . Table 8.) to calculate 
an initial T c g(ip) function from the input (V — Ic)(ip) curve. 
The actual accuracy of this formula is not an issue because a 
number of test runs show that the fitting process converges 
to the same solution T e ff(ip) curve when started from a wide 
range of initial functions i.e., the results are not sensitive 
to its choice. The solution T e s(<p) curve depends slightly on 
the choice of the initial V p (<p) function, but its effect is less 
than 10 — 30 K in any phase of the pulsation. 

2.2 The initial V p (tp) curve 

In contrast with the stability of the T e s(f) function, the situ- 
ation with the pulsational velocity, V p (ip), is quite different. 
We have found that, during the fitting process, the V p ((p) 
curve departs from a physically plausible shape if the fitting 
method is allowed to change its shape arbitrarily. Therefore, 
it is essential to find a reasonably good initial V p (<p) curve 
and, to get a reliable result, the solution V p (ip) curve has to 
be kept near to the initial one. 

2.2.1 The V p ((p) curve given by Liu (1991) 

iLiul (| 199 if ) published a template V p (ip) curve of RRab stars 
based on about a dozen observed Kad(<^) curves. The am- 
plitudes of the input curves were normalized to unity and 
they were phased to coincide at their minima. The average 
of these curves was the n cal c ulated to define the shape of 
a template V p (tp) curve. ILiul (|l99ll ) also derived a relation 
which gave the amplitude of the V p (fp) curve as a function 
of the V(ip) light curve amplitude with an accuracy of about 
5km/s. The actual template V p (ip) curve of an RRab vari- 
able is defined by the template shape and this amplitude for- 
mula. This template V p (ip) curve describes the input curves 
with about 5km/s accuracy, except on the steep descending 
branch between phases 0.85 and 1.0, where the accuracy is 
only about 25km/s. 

This average curve has, however, several drawbacks. 
Due to phase differences on the very steep descending 
branch, the template curve is l ess steep h ere as any of the 
individual curves (see Fig. 2 in [Liulll99ll ). For similar rea- 
son, the minimum is less sharp than most of the individual 
V p (ip) curves. 

Another problem with the template V p (tp) curve of lLiul 

|l99ll ) is its fixed shape. Our main purpose in developing 
the IP method was to study Blazhko stars showing large 



amplitude light curve modulation. Looking at the large va- 
riety of the light curves' shape of a Blazhko variable during 
the Blazhko cycle, it definitely does not seem to be a good 
decision to fix the shape of the V p (ip) curve to a unique tem- 
plate. Therefore, we have been looking for some connection 
between the V p (<p) curve, and the colour and light curves. 

2.2.2 V p (tp) curve calculated from the Ic(<p) light curve 

The only Kad observation that covers both the pulsation 
and modulation cycl e s of a Blazhko star was obtained by 
IChadid fc Chapellierl (|2006l ). Unfortunately, there is no si- 
multaneous spectroscopic and photometric observation of 
any Blazhko star that would permit to examine the possibil- 
ity of a light curve - V p curve connection valid for Blazhko 
stars in the course of the modulation. Therefore, in the fol- 
lowing we consider only unmodulated RRab stars. 

We used the pho t omet ric and V la ,a data of RRab stars 
collected by iKovac 1 (|2003h . From Kovacs's sample of 22 
RR Lyrae stars 20 are utilised, as RR Leo and W Crt has 
been omitted because of their discrepant curves. Addition- 
ally, we also used data of RV Oct. This sample contains all 
the RR Lyrae stars that have accurate Kad and photometric 
observations w ith good phase coverage of the pulsation. 

Following IKovac ] (|2003l ) we use p = 1.35 projection 
factor to transform radial velocities to pulsational velocities 
according to the formula: 

V p =p - V r ad, where p = 1.35. (5) 

Inspecting the observed K-ad(v) an d light and colour 
curves, we have found that the shape of the Ic(v) curve 
resembles the most to that of the V ra d(</>) curve. The V p vs. 
I c values of the 21 RRab stars are plotted in Fig.[T] Asterisk 
denotes that their phase- averaged value is subtracted from 
the data i.e., they are shifted to the same V p and magnitude 
zero points. This plot shows that the V p (<p) of an RRab star 
can be expressed as a function of its Iq(<p) brightness at any 
phase of the pulsation with sufficient accuracy. 

Statistical tests were performed to determine on how 
many parameters V p depends significantly and whether it 
is sufficient to express it with I c fc terms or it also depends 
on the colours. The r.m.s. of different degree regressions of 
V p * with terms of I c k , (V - Ic)', (B - V) 1 , I* c k ■ (V - Ic)', 
and I c fc • {B - V) 1 are plotted in Fig. [2] The plot reveals 
that no significant improvement can be achieved when we 
take into account a colour or compound term, too. It can 
also be seen in Fig. [2] that a 3rd order polynomial fit with 
Iq, Iq 2 , and Iq 3 terms describes V p with sufficient accuracy. 
However, if we fit the data directly by this polynomial, it will 
not be monotonic at the ends of the Iq range covered by 
the data. Therefore, to obtain a strictly monotonic V p - Iq 
relation, the following steps were performed. The plotted 
data were smoothed with a cubic spline and the spline was 
extrapolated linearly over the boundaries of the covered Iq 
range with about 0.05 mag at both ends. This spline was 
then fitted with a 7th degree polynomial. Finally, we found 
the following relation between V p and Iq, which represents 
the V p *(ip) curve of RRab stars similarly good or even better 
than Liu's template does. The formula: 

7 

V p *(v) = J2*- I c(v>) i (6) 

i=0 
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Table 1. Applied settings of the IP method. 



notation initial V p ((p) template V p (tp) weight factor 



Figure 1. V*(tp) vs. /£(</?) relation. The plotted values are cal- 
culated from the Fourier fits of the light and radial velocity curves 
of 21 stars sampled with 1/500 phase resolution. Asterisks denote 
that the time averages are subtracted from the data. The solid 
grey (cyan) line shows a 7th order polynomial fit according to 
Eq.© 




Figure 2. r.m.s. values of the regressions of V* with Iq 1 
I c ) 1 , (B - V) 1 , Jg fc ■ (V - IcY, and /£ fe • (B - V) 1 terms, 
circles denote the nth degree regressions when only Iq* 



(V- 
Filled 
terms 

are used. Open circles denote regressions where n — 1 terms of 
Iq° and one colour or compound term are taken into account. 
The r.m.s. does not drop significantly after the 3rd degree. The 
best 3rd order solution is achieved by a regression using the 
I'q , and 7q 3 terms. The arrow at the top left corner shows that 
data points here are at high r.m.s. values, between 9km/s and 
20km/s. 



where c .. 7 = (1.38, 194.8, -114.5, -1023, 

764,4220,-1660,-6700) 
fits the observations of the 21 stars with 3.51 km/s accuracy. 
The errors of the coefficients are: 

ct(co..t) = (0.08, 0.7, 3.7, 14, 40, 100, 110, 200) 



When running the IP code, initial V p (ip) curves derived 
either from Liu's template or from the V p - Iq relation are 
applied. As we are not sure that any of these V p (ip) curves 
are correct indeed, the shapes of the solution V p (<p) curves 
are not fixed to these templates. Unfortunately, as we men- 
tioned earlier the V p (ip) curve converges to a physically irre- 
alistic solution if it is allowed to change arbitrarily. Therefore 
the fitting algorithm treats the initial V p (tp) template as an- 
other curve that is fitted with an actual, varied V p (tp) curve. 
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large 
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Liu (19911 
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Liu (1991) 


large 







In this way, it becomes possible to give penalties to the tri- 
als depending on how much their V p (<p) curve deviates from 
the initial one. The rate of this penalty can be altered by a 
weight factor. Two weight factors are used. One allows V p (ip) 
to depart from the initial curve with about 4 km/s or less 
in any phase of the pulsation (medium weight). The other 
weight is 100 times larger, allowing practically no change 
in V p (ifi) (large weight). The variation of the initial V p ((p) 
curve is allowed with the addition of a 2nd order Fourier 
series, as our experiments show that higher orders lead to 
wavy solution curves. However, the results of the IP method 
depend more significantly on the quality of the input light 
curves than on the choice and variance of the initial V p (ip) 
curves. The 3 — 5 km/s accuracy of the initial Vp(<p) curves 
makes it possible to achieve similarly accurate results with 
the IP method using exculsively good quality photometric 
data as with earlier direct B-W analyses. 

As our method fits the output mathematical curves to 
input mathematical curves, the formal errors of the derived 
parameters according to the Levenberg-Marquardt fitting 
algorithm are meaningless and are independent from the 
uncertainties of the underlying observations. However, the 
results of the IP method are very sensitive to the accuracy 
and phase coverage of the photometric observations, since 
for good quality data the results depend only slightly on the 
internal settings of the code (Kad template, weight factor) 
while for worse data quality the results can show relatively 
large dispersion. Another problem of estimating the errors 
of the parameters derived from the IP method is that accu- 
rate fits of the input data can also be obtained for physically 
unreliable conditions, e.g, when the mass is not fixed by the 
pulsation equation, or when fitting the initial radial velocity 
curve has only a small weight. Therefore we regard the range 
of the output values from four different runs with different 
settings listed in Table 1 as an estimate of the inherent un- 
certainty of the IP method. 



3 TEST RESULTS 

The IP method has been tested on 9 RRab stars. We have 
selected test objects that have good BVIc light curves and 
Kad observations available. All of these objects were previ- 
ously subjected to B-W analyses by two or more authors. 
Four different runs of the IP method with different settings 
(detailed in Table [TJ have been performed for each test ob- 
ject to make uncertainity estimations possible. These set- 
tings are marked with the letters A, B, C, and D in the 're- 
mark/ref.' columns of Tables [2] and [4] Two different V p (tp) 
template curves (see Sects. 12.2.11 12.2.2|l have been applied 
with medium and large V p (tp) weight factors. 

The results of the tests are compared to the results of 
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Table 2. Comparison of the results of the IP method with literature data for the 9 test objects. The settings of the IP method using the 
notation of Table[T]are given in the 'remark/ref.' column. Italic numbers denote photometric metallicities and values that were calculated 
in this paper from the other parameters according to Eqs.[2][7] and [8] 



GCVS name (log P) 




[Fe/H] 


M vo 


d 


Ro 


TeffO 




log Sstat 


remark/ref. 






[mag] 


[PC] 


[Rq] 


[K] 


[an©] 




SW And (-0.354) 


A 




-0.10 


1.01 


504 


4.34 


6646 


0.56 


2.90 


13 




-0.10 


1.00 


507 


4.41 


6645 


0.57 


2.90 


c 




-0.10 


0.99 


508 


4.42 


6646 


0.57 


2.90 


D 




-0.10 


0.97 


513 


4.46 


6644 


0.58 


2.90 


McNamara & Feltz ( 


19771 


0.00 


0.90 


532 


4.45 


6680 


0.60 


2.92 


Cacciari et al. (1989 


) 


-0.15 


0.88 


537 


4.49 


6640 


0.60 


2.91 


Liu & Janes (19901 




-0.10 


0.97 


511 


4.36 


6517 


0.56 


2.91 


Jones et al. (19921 




-0.15 


1.10 


520 


4.06 








Fernlev (19941 




-0.15 


0.94 


522 










Kovacs (20031 




-0.10 


0.81 


555 


4.47 


6702 


0.58 


2.90 


Bono et al. (20031 




-0.24 


0.93 


523 










WY Ant (-0.241) 


A 




-1.55 


0.36 


1190 


6.38 


6527 


0.97 


2.82 


13 




-1.55 


0.34 


1200 


6.44 


6526 


0.99 


2.82 


C 




-1.55 


0.52 


1105 


5.93 


6525 


0.81 


2.80 


D 




-1.55 


0.49 


1119 


6.00 


6525 


0.83 


2.80 


Skillen et al. (19931 




-1.25 


0.63 


1035 


5.61 


6389 


0.72 


2.80 


Fernlev (19941 




-1.25 


0.55 


1089 










Kovacs (20031 




-1.55 


0.52 


1105 


6.10 


6296 


0.87 


2.81 


Bono et al. (20031 




-1.48 


0.48 


1124 










UU Cet (-0.217) 


A 




-1.10 


0.66 


1865 


5.72 


6420 


0.68 


2.75 


B 




-1.10 


0.62 


1896 


5.82 


6419 


0.71 


2.76 


G 




-1.10 


0.49 


2016 


6.19 


6420 


0.83 


2.78 


D 




-1.10 


0.44 


2061 


6.32 


6420 


0.88 


2.78 


Cacciari et al. (19921 IR Flux method 


-1.00 


0.63 


1887 


5.75 


6300 


0.70 


2.76 


Cacciari et al. (19921 SB method 


-1.00 


0.70 


1825 


5.36 




0.59 


2.75 


Cacciari et al. (1992 


) SB method 


-1.00 


0.52 


1982 


5.50 




0.62 


2.75 


Fernlev (19941 




-1.00 


0.62 


1899 










Kovacs (20031 




-1.10 


0.48 


2025 


6.25 


6258 


0.85 


2.77 


Bono et al. (20031 




-1.28 


0.56 


1949 










SU Dra (-0.180) 


A 




-1.60 


0.60 


695 


5.91 


6431 


0.63 


2.70 


B 




-1.60 


0.55 


712 


6.05 


6428 


0.69 


2.71 


C 




-1.60 


0.60 


696 


5.91 


6432 


0.65 


2.70 


D 




-1.60 


0.51 


724 


6.15 


6431 


0.71 


2.71 


A a 


[M/H] = -1.6; 


a/Fe] = 0.4 


0.57 


703 


5.96 


6435 


0.66 


2.71 


B tt 


[M/H] = -1.6; 


a/Fe] = 0.4 


0.51 


723 


6.13 


6433 


0.71 


2.71 


C a 


[M/H] = -1.6; 


a/Fe] = 0.4 


0.56 


707 


5.99 


6437 


0.67 


2.71 


D a 


[M/H] = -1.6; 


a/Fe] = 0.4 


0.47 


737 


6.24 


6438 


0.74 


2.72 


Liu & Janes (19901 




-1.60 


0.73 


640 


5.15 


6433 


0.47 


2.69 


Fernlev (19941 




-1.60 


0.63 


685 










Barcza (20031 




-1.60 


0.74 


647 


5.09 


6490 


0.66 


2.84 


Kovacs (20031 




-1.60 


0.62 


689 


5.85 


6293 


0.63 


2.70 


Bono et al. (20031 




-1.80 


0.26 


813 
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Table 3. Results (continued). 



GCVS name (log P) 


[Fc/H] 


M vo 


d 


Ro 


TeffO 


on 


logg s tat 


remark/ref. 




[mag] 


[PC] 


[«©] 


[K] 


P%] 




RR Gem (-0.401) 


A 


-0.15 


0.68 


1229 


4.74 


6882 


0.80 


2.99 


B 


-0.15 


0.66 


1243 


4.79 


6881 


0.82 


2.99 


C 


-0.15 


0.67 


1235 


4.76 


6884 


0.81 


2.99 


D 


-0.15 


0.65 


1247 


4.80 


6884 


0.83 


2.99 


Liu & Janes (1990) 


-0.20 


0.99 


1061 


4.05 


6699 


0.55 


2.96 


Fernlev (19941 


-0.30 


0.89 


1117 










Kovacs (20031 


-0.15 


0.74 


1196 


4.60 


6721 


0.74 


2.98 


Bono et al. (20031 


-0.29 


0.92 


1100 










TT Lyn (-0.224) 


A 


-1.20 


0.86 


626 


5.32 


6365 


0.58 


2.75 


B 


-1.20 


0.80 


643 


5.47 


6365 


0.62 


2.75 


C 


-1.20 


0.75 


660 


5.62 


6366 


0.66 


2.76 


D 


-1.20 


0.67 


685 


5.83 


6366 


0.73 


2.77 


Liu & Janes (19901 


-1.35 


0.75 


654 


5.40 


6284 


0.62 


2.77 


Fernlev (19941 


-1.35 


0.65 


690 










Kovacs (20031 


-1.20 


0.66 


686 


5.72 


6283 


0.69 


2.76 


Bono et al. (20031 


-1.56 


0.57 


717 










RV Oct (-0.243) 


A 


-1.15 


0.32 


1128 


6.35 


6567 


0.97 


2.82 


B 


-1.15 


0.28 


1149 


6.47 


6566 


1.02 


2.82 


C 


-1.15 


0.43 


1073 


6.04 


6565 


0.85 


2.81 


D 


-1.15 


0.40 


1089 


6.14 


6565 


0.89 


2.81 


A 


-1.75 


0.51 


1038 


5.94 


6538 


0.82 


2.80 


B 


-1.75 


0.43 


1076 


6.16 


6537 


0.90 


2.81 


C 


-1.75 


0.61 


990 


5.67 


6537 


0.73 


2.79 


D 


-1.75 


0.56 


1014 


5.81 


6536 


0.77 


2.80 


Skillen et al. (19931 


-1.75 


0.76 


905 


5.26 


6437 


0.62 


2.79 


Fernlev (19941 


-1.75 


0.68 


956 










Bono et al. (20031 


-1.71 


0.37 


1102 










AR Per (-0.371) 


A 


-0.05 


1.02 


489 


4.25 


6719 


0.54 


2.92 


B 


-0.05 


1.03 


487 


4.23 


6720 


0.54 


2.91 


C 


-0.05 


1.11 


470 


4.09 


6715 


0.49 


2.91 


D 


-0.05 


1.11 


469 


4.08 


6716 


0.49 


2.91 


Liu & Janes (19901 


-0.30 


0.97 


500 


4.15 


6672 


0.53 


2.93 


Fernlev (19941 


-0.30 


0.87 


524 










Kovacs (20031 


-0.05 


0.72 


561 


4.45 


6833 


0.61 


2.93 


Bono et al. (20031 


-0.30 


0.92 


512 










BB Pup (-0.318) 


A 


-0.45 


0.80 


1648 


4.65 


6790 


0.57 


2.86 


B 


-O.45 


0.76 


1677 


4.73 


6789 


0.60 


2.86 


C 


-0.45 


0.82 


1633 


4.61 


6793 


0.55 


2.85 


D 


-0.45 


0.79 


1658 


4.67 


6794 


0.58 


2.86 


Skillen et al. (19931 


-0.60 


1.21 


1329 


3.97 


6551 


0.40 


2.84 


Fernlev (19941 


-0.60 


1.13 


1416 










Kovacs (20031 


-0.45 


0.87 


1596 


4.40 


6712 


0.49 


2.84 


Bono et al. (20031 


-0.64 


0.65 


1766 
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Figure 3. Derived light curves (solid lines), input light curves (dashed line) and actual photometric measurements (dots) of RR Gem, 
SW And and UU Cet, stars representing two better and a poorer quality input light curves, respectively. The four runs of the IP method 
with different settings are shown with different colours, however, the differences between these solutions are in most cases less than the 
thickness of the lines. 
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Figure 4. Differences between the initial and solution light and 
colour curves for RR Gem. Results of the IP method with four 
different settings are plotted with different colours. The periodic 
nature of the residuals arises from differences of the order of the 
Fourier sums used to describe the input and the fitted output 
curves. 



earlier B-W analyses and to the Mvo absolut e brightness 
value s derived from pulsational luminosities (|Bono et al.l 
I2003T I in Table [2] This table lists the source of the data, 
the metallicity, mean absolute visual magnitude (Mvo), dis- 
tance (d), mean radius (Ro), mean effective surface temper- 
ature (Tefio), mass (9Jt), and static gravity (log g s t a t) of the 
stars. If the mass, logpstat, and/or distance were not given 
in the original paper cited, they were calculated according 
to Eqs. El and [S] respectively. These values are typeset 
with italics in Table [2] The static gravity and distance are 
expressed as: 




0.4 0.6 
Phase 

Figure 6. Radius variation curves of SW And derived from the IP 
method compared to results published by other authors. Continu- 
ous lines denote the solutions of the IP method with four different 
settings. Symbols denote literature data: triangles - Cacciari et 
al. (1989), filled circles - Liu &; Janes (1990), open circles - Jones 
et al. (1992). 



log fifstat = log M - 2 log R + 4.438. 
and 

logd = 0.2 ■ (mvo - Mvo - Av + 5) 



(7) 



(8) 



where mvo is the mean apparent V magnitude and Av is 
the interstellar extinction in V band. 

The original B-W analyses used spectroscopic metallic- 
ity values except iKovacj (|2003l ) who calculated the [Fe/H] of 
the studied variables from the light curve Fourier parameters 
(|jurcsik &c Kova cs 1996, Eq. 3, photometric metallicity). In 
Table [2] [Fe/H] values corresponding to this formula are type- 
set with italics. As the photometric metallicities were cali- 
brated on spectroscopic data the two metallicities are iden- 
tical within the limits of the uncertainties. For all stars, with 
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Figure 5. Solutions for the V ra d(¥>) curves of the 9 test objects are compared to their V r * d measurements. Gray (green) lines show Liu's 
template, black (red) lines mark the Vrad(v) template according to Eq. \E\ The medium weight solutions are drawn with solid lines, the 
large weight solutions are plotted with dashed lines. Dots show the actual V* ad measurements with mean values subtracted. Error bars in 
the bottom left corners of each panel indicate the mean errors of the observations. V r * d curves derived from Iq light curves are somewhat 
wavy when the light curve is of poorer quality and its high enough order Fourier fit is not smooth. Our V p (ip) curves are transformed to 
V ra( j(ip) using the actual value of the projection factor (p) given in the original papers. The centre-of-mass velocities of the observations 
are subtracted from the data. The slight differences between the observed radial velocity curves and the solution pulsation velocity curves 
tra nsformed to radial ve lociti es may either or i ginat e from the pulsation phase and period dependence of the transformation factor p (see 
e.g. ISabbev et ah! l ll995l) and lNardetto et al.l ll2007t) for Cepheids) or from minor defects of the V ra d templates used. 



the exception of RV Oct, the spectroscopic and photometric 
metallicities agree within 0.0 — 0.3 dex. The spectroscopic 
metallicity of RV Oct is —1.75 dex, while photometrically 
— 1.15 dex is derived. Therefore, we performed all the 4 runs 
of the method with both metallicity values for this star. For 
other test objects we accepte d the photometric metallici- 
ties. The metallicities u sed by Bono et al. ([20031 ) are taken 
from the compilation bv lFernlev et al.f "( 1998T ). We have also 



checked how the results change if a enhanced models are 
used. The IP method was run for SU Dra, for one of the 
metal poor stars using [M/H] = —1.6, [o/Fe] = 0.4 atmo- 
sphere models, too. In Table [2] these solutions are denoted 
by 'a'. The effect of using a enhanced models are small, it 
is about +5K and —0.04 mag in T £f and My, respectively. 
For compatibility reasons (all the previous direct B-W anal- 
yses used solar scaled atmosphere models) in the following, 
solutions of the IP method using solar scaled atmosphere 
models are regarded. 

The derived light curves reproduce the input light 
curves to an accuracy of 0.01 — 0.03 mag which is about 
the observational uncertainty as shown in Figs [3] and [4] The 



results of the fitting process are shown by the examples of 
SW And and UU Cet in Fig. [3] The input light curves of 
SW And are well defined by the observations, especially in 
V band, while UU Cet has poorer photometric data, partic- 
ularly at the lower part of the rising branch. The method 
produces more stable result (i.e., the difference between the 
four solutions are smaller) for those stars that have better 
quality light curves. Fig. [4] shows the difference between the 
input and derived light and colour curves of RR Gem. The 
differences are less than about 0.01 mag in V and (V — 7c) 
and less than 0.02 mag in (B — V). 

As a result of the IP method, solution V p (ip) curves 
are also derived, which can be converted to Kad(</?) curves. 
These curves are compared to the actual Kad measurements 
of our 9 test objects in Fig. [5] Our V p (ip) curves were trans- 
formed to Urad(¥>) using the p transformation factors given 
by authors of the measurements. The radial velocity data 
have typically 1-4 km/s uncertaintie s, as error bar s indi- 
cate in the plots. Data are taken fro m ICacciari et al.l(ll987l. 
SW A nd), Iciementini et al.l (|l99fj| , UU Cet). I Liu fc Janes 
(| 19891 . SU Dra, RR Gem, TT Lyn, AR Per) and lSkillen et al 
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-0.4 -0.35 -0.3 -0.25 -0.2 -1.5 -1 -0.5 

log P [Fe/H] 



Figure 7. log P vs. log g s tat values from Table [2] The results of 
the IP method are plotted with black dots while literature data 
are marked with gray colour. The scatter of the literature data 
and the data from the four runs of the IP method with different 
settings indicate the inherent uncertainties of direct B-W analyses 
and of the IP m ethod. Solid l ine corresponds to the empirical 
formula given bv lJurcsikl lll998l . Eq. 15). 



Figure 8. [Fe/H] vs. Myo values from Table [2] Our results are 
marked with black symbols while literature data are plotted with 
gray colour using the same symbols. Solid line shows the zero age 
horizontal branch relation given bv lSandage &c T ammann (|2006l . 
Eq. 9). 



(1993, WY Ant, RV Oct, BB Pup). Note that the Kad ob- 
servations were not taken into account by the IP method 
in any way. In most cases the radial velocity curves derived 
from the IP method agree within the uncertainties with the 
observations. There are some systematic differences within 
about 5km/s range apparent, most probaly due to minor 
defects of both Liu's V ra A template and the V p — Ic ralation. 
The largest deviations are observed at around maximum 
brightness phases of UU Cet. However, the good agreement 
between the physical parameters derived even for this star 
from the IP method and from direct B-W analyses shows 
that our method is not sensitive to these small inaccuracies 
of the solutions for the radial velocity curves. 

From the solution V p (<p) curve and from the mean ra- 
dius the variation of the radius with phase, R(f), can be 
derived. The four solution R(<fi) curves of SW And are com- 
pared to the curves given by authors of earlier B-W analyses 
in Fig. [6] The agreement here is fully satisfactory. 

The physical parameters derived from the IP method 
are in reasonably good agreement with results of other meth- 
ods, as shown in Table[2] For individual stars the Mv o values 
derived by different authors differ by about 0.1 — 0.3 mag, 
while the different settings of the IP method give results 
within 0.2 mag range. As there is no objective reason to de- 
cide which of the literature data are the most correct, the 
direct comparison of our results with them is not obvious. 
E.g., the absolute magnitudes derived from the IP method 
for B B Pup agree well w i th th e val ue given by Kovacsl 
(|2003l ). while ISkillen et all (|l993l) and iFernlevI <|l994) give 
0.2 — 0.3 mag fainter results. 

Therefore, we decided to check the consistency of our 
results with empirical and theoretical relations valid for 
RR Lyrae variables in order to become certain about the 
reliability of our results. 

There is a well-defined relation betw een the log P and 
log g values of pulsating variables (see e.g. lFern"ielll995l ). The 
log P vs. log (?atat points of the IP method solutions and the 
data given by other authors for the test objects are plot- 




0.2 0.4 
0.471 -1.132 log P + 0.190 log Z 

Figure 9. Absolute Iq brightness values (MI0) derived from 
the IP method (black symbols) and from direct B-W analy- 
ses (grey symbols) compared to the ir predicted values f rom the 
M I0 (logP,logZ) relation denned bv lCatelan et al.l d20u4 Eq. 3). 
Solid line shows the x = y equality. 



te d in Fig, ffj Sol id line shows the empirical relation given 
by Ijurcsikl (|l998l . Eq. 15) . Our results fit to the empirical 
relation similarly well as the data of direct B-W analyses. 

Another important relation exists between the metallic- 
ity and the average absolute visual brightness of Horizontal 
Branch stars. The [Fe/H] vs. Myo results of the IP method 
and of other authors are compared to zero age horizontal 
bran ch (ZAHB) model predictions (|Sandage fc Tammannl 
120061 ) in Fig. [8] This fi gure also shows that our results are in 
good accordance with theoretical predictions. The ambigu- 
ity of our results arising from the different settings of the IP 
method is smaller for most of the stars than the ambiguity 
between the results of earlier analyses. 

The reliability of the results was also checked and com- 
pared with direct B-W data using PL and PLC relations of 
the absolute Ic brightness as shown in Figs. and 1101 The 
Mi (logP,logZ) and the log P (I C , V - Ic) relati o ns de fined 
by ICatelan et all (|2004) and iKovacs fc Walker! (|200ll ) are 
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Figure 10. Period — V brightness— (V — Ic) colour realtion ac- 
cording to Eq. 8 of iKovacs fc Walker! 1 120011) . The constant of 
this equation (C) was calculated by least squares fits to our data 
(black dots) and for the literature data (grey dots). Within the 
uncertainties the same constant of the PhC(V, 7c) relation was 
obtained from both data sets with the same r.m.s. scatter of 
0.03 mag. Solid and dashed lines correspond to the y = x + C 
equalities for IP and literature B-W data, respectively. 

used for these tests. Again, no discrepancy of the results of 
the IP method are evident. 



4 APPLICABILITY 

We have also tested how the IP method works if not all the 
necessary input data are available. The results are shown in 
Table [4] for SW And and RR Gem, the two stars that have 
the best quality photometric observations among the test 
objects. 

4.1 Two colour (Vic or BV) light curves 

The IP method was tested using only V and Ic or B and 
V light curves. As the initial T c r (ip) and V v (ip) curves were 
defined in Sect. 2 partially using the 7c data, if only BV 
data are available the initial T e ff(tp) curve is defined from 
empirical T e ff(73 - V, [Fe/H]) relation (Table 2., Eq. 3 in 
lAlonso et al.lll999l ). and in this case, Liu's template V p [ip) 
curve can only be used. 

The results of four and two runs with different settings 
taking only V and 7c or B and V light curves of SW And 
and RR Gem into account are given in Table [4] marked with 
'vi' and 'bv' in the 'remark' column. 

The IP method using only V and 7c light curves gives 
0.1 mag brighter and fainter Mvo values, 0.2 7?© larger and 
smaller mean radii, and 0.1 9JIq greater and smaller masses 
than the BVIc solutions for SW And and RR Gem, re- 
spectively, while the temperature of both stars are cooler by 
20 — 40 K according to these solutions. These values are not 
significantly discrepant from literature data, the agreement 
for RR Gem is even better than when using BVIc data. 

If only B and V data are utilized the temperature so- 
lution of both stars are hotter by 30 — 40 K. In this case the 
other parameters of SW And, though are still resonable, dif- 
fer from the direct B-W results more than the results of the 



BVIc solutions. The BV solutions give, however, unreliably 
bright Mvo and large radius and mass for RR Gem. 

Therfore, we conclude that the IP method give accept- 
able solutions only if at least V and 7c band light curves are 
available, but for B and V data it may lead to unreliable 
results. In direct B-W analyses it was also found by many 
authors that light curves observed in longer wavelength pass- 
bands (VIcJHK) give more accurate results than those ob- 
served in bluer passbands (UB). 

4.2 Unknown zero points of the colour curves 

To fix the zero points of the light and colour curves the red- 
dening and standard magnitudes of the variable and/or the 
comparison stars have to be known. However, these condi- 
tions are not always met. 

We have also performed test runs to check the results 
of the IP method when the (73 — V)o and (V — Ic)o zero 
points are not known. During these tests the colour curves' 
zero points were not fixed but they were also fitted by the 
algorithm. Our tests show that, for good quality light curves, 
the algorithm converges to the appropriate zero points. Even 
when the initial values of the zero points of the colour curves 
are off about 0.1 mag from their correct values, the method 
finds them within about 0.01 — 0.03 mag accuracy. These 
differences cause only marginal changes in Mvo, R, an d , 9JI, 
0.06 mag, 0.06 Rq, and 0.02 2ft©, respectively. The output of 
these four runs are listed in Table[4] marked with 'zp' in the 
'remark' column. For these results the zero point differences 
are also given in the last two columns with respect to the 
observed zero points (A(B — V)o and A(V — /c)o)- 

Finally, we have also checked what results are obtained 
if neither B observations, nor (V — 7c)o zero point were 
available. These results are denoted by 'vi zp' in Table [4] 
For the two stars even these solutions remain reliable. We 
have to emphasize, however, that in the case of incomplete 
input data the quality of the observed light curves becomes 
even more important. For poorer quality light curves these 
limited solutions may lead to physically meaningless results. 



5 SUMMARY AND CONCLUSIONS 

We introduced a new Inverse Photometric Baade-Wesselink 
method to determine physical parameters of RRab stars ex- 
clusively from their multicolour light curves. It was shown 
through the example of 9 test objects that for good quality 
BVIc light curves the IP method gives similarly accurate 
results to other analyses, without the need for spectroscopic 
K-ad observations. 

The success of the IP method depends strongly on the 
quality of the light curves. For the best quality input curves 
the method is applicable even when only V and 7c observa- 
tions are available and when the zero points of the colours 
are not known. Note that any error in mvo influences only 
the distance determination. 

The method is an inverse method, since the primary 
curves are the V p (<p) and T e g(ip) functions and the output 
light curves are derived from them through the application of 
static atmosphere models, while direct B-W methods derive 
Teff(vs) and \ogg c ff(tp) curves directly from the observations. 
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Table 4. Results of the applicability tests on SW And and RR Gem. The two-colour (Vic and BV) solutions are marked with 'vi' and 
'bv', respectively, and solutions when the zero points of the colours are also fitted are marked with 'zp' in the 'remark' column. 



GCVS name (log P) 


[Fe/H] 


M vo 


d 


Ro 


T e ff0 


m 


loggstat 


A(B — V)o 


A(V - J)o 


remark* 




[mag] 


[pc] 


[Re] 


[K] 






[mag] 


[mag] 


SW And (-0.354) 


A vi 


-0.10 


0.92 


527 


4.63 


6609 


0.64 


2.91 






B vi 


-0.10 


0.91 


529 


4.65 


6609 


0.65 


2.91 






C vi 


-0.10 


0.89 


534 


4.69 


6609 


0.66 


2.91 






D vi 


-0.10 


0.91 


531 


4.66 


6608 


0.65 


2.91 






C bv 


-0.10 


1.14 


475 


4.08 


6681 


0.46 


2.88 






D bv 


-0.10 


1.13 


479 


4.11 


6680 


0.47 


2.88 






A bvi zp 


-0.10 


0.95 


520 


4.42 


6707 


0.57 


2.90 


0.005 


0.027 


B bvi zp 


-0.10 


0.98 


512 


4.41 


6665 


0.56 


2.90 


-0.005 


0.016 


C bvi zp 


-0.10 


0.93 


525 


4.47 


6705 


0.58 


2.90 


0.005 


0.026 


D bvi zp 


-0.10 


0.92 


527 


4.48 


6708 


0.59 


2.90 


0.006 


0.027 


A vi zp 


-0.10 


0.69 


584 


4.75 


6842 


0.68 


2.92 




0.062 


B vi zp 


-0.10 


0.70 


583 


4.75 


6841 


0.68 


2.92 




0.061 


C vi zp 


-0.10 


0.68 


589 


4.80 


6840 


0.70 


2.92 




0.061 


D vi zp 


-0.10 


0.67 


591 


4.81 


6842 


0.70 


2.92 




0.062 


RR Gem (-0.401) 


A vi 


-0.15 


0.80 


1168 


4.52 


6858 


0.71 


2.98 






B vi 


-0.15 


0.79 


1171 


4.54 


6856 


0.71 


2.98 






C vi 


-0.15 


0.77 


1183 


4.58 


6859 


0.73 


2.98 






D vi 


-0.15 


0.77 


1183 


4.58 


6859 


0.73 


2.98 






C bv 


-0.15 


0.40 


1400 


5.34 


6914 


1.08 


3.02 






D bv 


-0.15 


0.40 


1401 


5.34 


6914 


1.08 


3.02 






A bvi zp 


-0.15 


0.66 


1241 


4.74 


6907 


0.80 


2.99 


0.000 


0.012 


B bvi zp 


-0.15 


0.63 


1258 


4.83 


6890 


0.84 


2.99 


-0.004 


0.008 


C bvi zp 


-0.15 


0.63 


1259 


4.76 


6940 


0.81 


2.99 


0.007 


0.019 


D bvi zp 


-0.15 


0.64 


1255 


4.81 


6893 


0.83 


2.99 


-0.004 


0.008 


A vi zp 


-0.15 


0.74 


1199 


4.52 


6946 


0.71 


2.98 




0.022 


B vi zp 


-0.15 


0.72 


1205 


4.51 


6973 


0.70 


2.98 




0.029 


C vi zp 


-0.15 


0.73 


1202 


4.57 


6918 


0.73 


2.98 




0.015 


D vi zp 


-0.15 


0.77 


1183 


4.57 


6865 


0.73 


2.98 




0.002 



*: On the notation A, B, C, and D see Table [T] 



The IP method varies these primary curves during a nonlin- 
ear least squares fitting process so that the derived light and 
colour curves best match the observed input curves. The ini- 
tial T e g(ip) curve is calculated fro m the (V — Ic)(<p) colour 
curve using an empirical formula (jBessel et al ] 1 19981 . Table 
8.). The initial V p (ip) curve is defined twofold. We have de- 
rived a new 7c - V p relation which yields the V p (<p) curve 
from the Ic((p) light curve with an accurac y of 3.5 km /s 
(Eq. 0. Alternatively, the V p (ip) template of iLiul (|l99ll ) is 
used. 

The method was developed to extract physical informa- 
tion from accurate multicolour RR Lyrae light curves. Dur- 
ing the development of the IP method special attention was 
payed to its applicability to RR Lyrae stars showing Blazhko 
modulation. The applicability and the results of the method 
fo r a strongly modulate d RR Lyrae star, MW Lyr, is shown 
in ljurcsik et al] (|2008bh . The IP method is the only possibil- 
ity today to derive the variations of the physical parameters 
of a Blazhko star throughout the modulation cycle. 

Another field of interest of the application of this new 



method is to study RRab stars in globular clusters or dwarf 
galaxies since there are very few radial velocity observations 
made of these faint variables but, thanks to the CCD tech- 
nique, accurate multicolour data are available. 

The IP method has been specialised to fundamental 
mode RR Lyrae stars, although with the appropriate choice 
of the initial V p (ip) curve and using adequate pulsation equa- 
tion it could be adapted for other type of radial pulsators. 
Using appropriate model atmosphere tables, the method is 
applicable for other passbands, as well. To do this, how- 
ever, empirical formulae between the appropriate colours 
and brightnesses, and the V p (ip) and T e g(ip) curves have to 
be defined first. 
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